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ABSTRACT 

Constraints on the ionization structure of the Intergalactic Medium are derived as di- 
rectly imposed by observations in conjunction with the results of numerical simulations 
for structure formation. Under the assumption that the population of sources dominat- 
ing the UV ionizing background at z < 6 is the same population which reionized the 
IGM, it is shown that consistency with measurements of the mean Lya transmitted 
flux at high redshifts requires the epoch of hydrogen reionization of the IGM to have 
occurred at a redshift z T \ < 11, independent of the space density of the sources. The 
upper limit on the reionization redshift depends only on the shape of the UV spectra of 
the sources. Consistency with constraints on the reionization epoch from the Wilkinson 
Microwave Anisotropy Probe requires the sources of photoionization to have had hard 
spectra, such as QSOs or Population III stars. The only way to escape these conclu- 
sions is either: 1. the sources which dominate the photoionization background at z < 6 
are the remnants of a population that was a much more prodigious source of ionizing 
photons at earlier times, or 2. the sources responsible for the photoionization of the 
IGM are an unknown population that contributes negligibly to the ionization of the 
IGM at z < 6. The evolution of the QSO luminosity function over the range 3 < z < 6 
is estimated from recent QSO counts, and the fraction of ionizing photons arising in 
QSOs is evaluated under the assumptions of either pure luminosity evolution or pure 
density evolution. It is shown that QSOs dominate the UV ionizing background for 
z < 3.5, but that it is unlikely that more than one half to one third of the ionizing 
background radiation originates in QSOs in the redshift range 4.5 < z < 6. QSOs 
acting alone could not have reionized the hydrogen in the IGM prior to z « 4. If the 
QSOs had hard spectra, however, they may have reionized the helium in the IGM as 
early as z w 5. The possibility that the IGM was reionized by low luminosity AGN is 
discussed. 
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1 INTRODUCTION 

The ionization structure of the Intergalactic Medium (IGM) 
at high redshift has been the focus of much recent atten- 
tion. On the basis of measurements of the mean transmitted 
Lya flux in spectra of the highest known redshift Quasi- 
Stellar Objects (QSOs) by the Sloan Digital Sky Survey 
(SDSS), it has been suggested that the reionization of the 
IGM had just completed at z > 6 (Becker et al. 2001). By 
contrast, measurements of fluctuations in the Cosmic Mi- 
crowave Background (CMB) by the Wilkinson Microwave 
Anisotropy Probe ( WMAP) suggest the IGM was reionized 
earlier, at z > 11 (2a lower limit) (Kogut et al. 2003). The 
most obvious candidate sources of reionization are stars and 
QSOs. An assessment based on the results of the Hubble 
Ultra-Deep Field (UDF) and the Great Observatories Ori- 



gins Deep Survey (GOODS) by Bunker et al. (2004) sug- 
gests that the galaxies identified at z ~ 6 are too few to ac- 
count for the reionization of the IGM by z = 6. On the other 
hand, Yan & Windhorst (2004) and Stiavelli, Fall & Panagia 
(2004) dissent from this view, arguing that, with plausible 
assumptions, these sources are adequate. Yan & Windhorst 
further argue that galaxies appear to be the only sources of 
reionization. They claim QSOs contribute negligibly to the 
overall ionizing photon budget, falling short of the required 
numbers by more than two orders of magnitude. Meiksin 
& White (2004) argue instead for a much smaller discrep- 
ancy for the quasars. The range of views in the literature 
highlights the sensitivity of conclusions regarding reioniza- 
tion to prevailing uncertainties, in particular the luminosity 
function of the sources, their nature (especially their emer- 
gent UV spectra), and the structure of the IGM. 
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A question distinct from the sources responsible for 
reionization is the nature of the sources which dominate the 
UV photoionizing background at z < 6, when the IGM is 
clearly ionized. The populations need not be the same, as the 
ionizing contribution from the sources responsible for reion- 
izing the IGM may decline with time, while other sources be- 
gin to dominate the UV background of the post-reionization 
universe. For instance, the mean escape fraction of ionizing 
photons from galaxies at high redshifts is currently unclear, 
and could in principle be sufficiently high to enable galaxies 
to reionize the IGM. Steidel et al. (2001) report a detection 
of UV ionizing photons from the the brightest Lyman break 
galaxies in a sample at z ~ 3, suggesting an escape fraction 
of about 10%, although this difficult measurement has not 
yet received independent confirmation. Attempts to quantify 
the mean escape fraction have generally resulted only in up- 
per limits. For instance, Fernandez-Soto, Lanzetta & Chen 
(2002) constrain the mean escape fraction for galaxies in the 
range 1.9 < z < 3.5 to less that 4%, which may be too small 
for galaxies to maintain the ionization of the IGM. Applying 
an escape fraction of ~ 10% to the full range of galaxy lumi- 
nosities, Steidel et al. obtain a value for the global emissivity 
of ionizing radiation that is inconsistent with the emissivity 
required by numerical simulations of the IGM to reproduce 
the measured values of the mean transmitted Lya flux at 
2Ri3. The high emissivity obtained by Steidel et al. would 
result in a much higher Lya transmitted flux through the 
IGM than measured, unless the sources were very short-lived 
(eg, bursting at just z ^ 3) (Meiksin & White 2003). This 
suggests that the escape fraction found by Steidel et al. may 
be anomalously high. While an escape fraction of 1 — 3% may 
be just sufficient for galaxies to maintain the ionization of 
the IGM at z ~ 3, it is unclear whether this fraction is suf- 
ficient for galaxies to meet the requirements for reionizing 
the IGM (see Stiavalli et al. ). 

Constraints on the epoch of reionization and on the level 
of ionization in the post-reionization IGM have usually been 
treated as separate issues. In this paper, it is shown that the 
two are closely related. Any model that meets the require- 
ments for the reionization epoch implies a level of ioniza- 
tion of the IGM just after reionization is complete. This 
level must be consistent with constraints imposed by mea- 
surements of the mean Lya transmitted flux through the 
IGM. In particular, if the high Lya optical depth measure- 
ments at z « 6 designate the tail end of the reionization 
epoch, then the rate of UV ionizing photon production by 
the sources must 'gracefully exit' the reionization epoch at 
the level required to match the measured Lya optical depth 
measurements at z 6. 

A second purpose of this paper is to provide a statistical 
assessment of the contribution of QSOs to the photoioniz- 
ing background at z > 3, and to the overall budget of ion- 
izing photons at z ~ 6, using the current best estimates of 
the faint (Hunt et al. 2004) and bright (Fan et al. 2001, 
2004) ends of the QSO luminosity- redshift distribution. Es- 
timates of the UV background based on QSO counts arc 
not very well constrained at high redshifts (z > 3) due to 
the low sampling of the luminosity-redshift plane by existing 
surveys. The uncertainty is sufficiently large to allow for a 
substantial contribution to the UV metagalactic background 
from QSOs at high redshifts (Meiksin & White 2004). 



2 THE IONIZATION POROSITY OF THE IGM 

The epoch of reionization by sources of number density ns 
may be characterised by the time at which the porosity 
Q — ns{Vi) exceeds unity, where (Vi) is the average volume 
of the universe ionized per source. At earlier stages, the fill- 
ing factor of the remaining neutral gas is / = 1 — Q. For 
any given evolutionary scenario of the sources, the porosity 
is determined by the total rate of production Ns of ionising 
photons per unit volume, which in turn scales like the total 
emissivity of the sources. On the other hand, for a uniform 
medium, the neutral fraction, and hence the Lya optical 
depth r a , of the IGM scales inversely with the total source 
emissivity. As a consequence, the product of the two will re- 
flect only the evolution and spectral shape of the sources 
(and the properties of the IGM), but be independent of 
their actual numbers. For any given source history and spec- 
tral shape, the Lya optical depth just following the epoch 
of reionization is thus well-specified and must be consis- 
tent with measured values. Turning this argument around, 
if the evolution of the emissivity and attenuation length 
are known, or reasonable guesses made based on observa- 
tions or numerical simulations, then measurements of the 
Lya optical depth may be used to provide a direct estimate 
of the porosity of the IGM, and hence of the epoch of reion- 
ization. 

In reality, the IGM is not uniform, but clumped into 
the Lya forest, so that the relation between the porosity and 
the measured mean Lya optical depth is indirect. (Here, the 
mean Lya optical depth is defined as f a = — log(/ a ), where 
(fa) is the mean transmitted Lya flux from a background 
source.) It is necessary to appeal to models of structure for- 
mation to establish the relation. Within the context of the 
models, the intimate relation between the porosity and mean 
transmitted Lya flux persists. Setting the source emissivity 
at a level to reionize the IGM at a given epoch fixes the 
post-reionization photoionization rate of the IGM, and so 
determines the mean transmitted Lya flux just after reion- 
ization is completed (for any given cosmological model) . Any 
reionization model must yield predictions consistent with 
these measured values. For instance, if the predicted pho- 
toionization rate were too high, the predicted mean trans- 
mitted Lya flux levels would exceed those measured. Con- 
straints on the reionization epoch and the post-reionization 
Lya transmitted flux values must be satisfied simultane- 
ously. 

The results of a ACDM simulation run to investigate 
effects of radiative transfer on the Lya forest (Meiksin & 
White 2004) are used to quantify the emissivity required to 
maintain the ionization of the IGM over the redshift range 

3 < z < 6 at a level consistent with current measurements of 
the mean transmitted Lya flux. The simulation was run us- 
ing a pure Particle Mesh (PM) dark matter code, and it was 
assumed the gas and dark matter have the same spatial dis- 
tribution, which is found to give an adequate description of 
the structure of the IGM (Meiksin & White 2001, 2003). 
The parameters used for the simulation are Q.m = 0.30, 
n A = 0.70, n b = 0.045, h = fVlOOkms^Mpc- 1 ^ 0.70, 
as = 0.92 and slope of the primordial density perturbation 
power spectrum n = 0.95. These values will be used through- 
out. The model is consistent with existing large-scale struc- 
ture, Lya forest flux distribution, cluster abundance and 



© 0000 RAS, MNRAS 000, 000-000 



Constraints on the ionization sources of the high redshift intergalactic medium 3 



a 6 — 




Figure 1. Mctagalactic emissivity (proper) at the Lyman edge 
(in units of 10 26 ergs -1 Hz -1 Mpc~ 3 ), as a function of redshift. 
The data points (filled circles and upper limit point) are the emis- 
sivity values required in a ACDM model of the IGM to match the 
measured values of the transmitted Lya flux through the IGM. 
The shaded regions show the la (dark) and 2a (light) ranges for 
a power-law fit (solid line) to the evolution in emissivity. The long 
dashed line shows the required source emissivity, allowing for a 
boost in the total emissivity by 50% from diffuse emission by the 
IGM. The squares show the estimated contribution from QSOs 
assuming either a pure luminosity evolution model (filled squares) 
or a pure density evolution model (open squares) for the QSO lu- 
minosity function. (See § 3.) The short dashed and dotted lines 
show the deficit left by these (respective) QSO models from the 
required source emissivity (long dashed line). The deficit must 
be made up by sources other than the QSOs. (A mctagalactic 
spectral index of qmg = 1 nas been assumed.) 



WMAP constraints (Meiksin & White 2004). The simula- 
tion was run using 512 3 particles and a 1024 3 force mesh, 
in a cubic box with (comoving) side length 25/i~ 1 Mpc, ad- 
equate for obtaining converged estimates of the Lya pixel 
flux distribution, the photon attenuation length and the re- 
quired mean emissivity at the Lyman edge (Meiksin & White 
2003, 2004). A description of the parallel PM code is given 
in Meiksin & White (2003). 

In the limit of a medium with a photon attenuation 
length at the Lyman edge short compared with the horizon 
scale, the photoionization rate Y is related to the emissivity 
eL at the Lyman edge by 

L ~ , rr (1) 

hp{3 + O.MG) 

where ol is the photoeletric cross-section at the threshold 
frequency ul, hp is the Planck constant, ro is the photon 
attenuation length, and the spectrum of the ionizing back- 
ground is taken to vary at frequencies above the threshold 
as a power-law e v = /vL,)~ aMG ■ In Fig. 1, the resulting 
constraints on the metagalactic emissivity using the data 
in Table 3 of Meiksin & White (2004) are shown. These 



data were based on ACDM model requirements to match 
the mean transmitted Lya flux measurements from a vari- 
ety of sources (see Meiksin & White 2004, Appendix B). It is 
emphasized that the emissivity and attenuation length are 
solved for self-consistently in the model. Simulation-based 
estimates of the required emissivity in the literature usu- 
ally rely on an assumed attenuation length (or, equivalenty, 
effective optical depth at the Lyman edge) based on sta- 
tistical representations of the measured Lya forest absorber 
properties (like H i column density and Doppler parameter) . 
Uncertainties in the measured distributions, however, pro- 
duce uncertainties in the derived attenuation length as large 
as a factor of 2-3 (Meiksin & Madau 1993, Meiksin & White 
2003, 2004). This uncertainty translates directly into a com- 
parable uncertainty in the required source emissivity. For 
this reason, it is crucial that the inferred emissivity from a 
simulation be consistent with the attenuation length that 
the same simulation predicts. 

A least-squares fit to the emissivity values from Meiksin 
& White (2004) over 2.75 < z < 5.5 gives the comoving 
emissivity 



ul m Ayio + ^ MG ) (1 + z)"' feergss 1 Hz 1 Mpc 
where ,4mg = x 10 25 and 7 = - 



(2) 



-1.6 



0.61og(A M G/8.4 x 10 25 ). The errors on Amg are the la and 
2<r limits. (The fit is acceptable at the 13% confidence level 
according to the y^-test.) The upper limit at z « 6 is con- 
sistent with this fit, as shown in Fig. 1. 

The total rate at which ionizing photons are generated 
by sources, per unit comoving volume, is related to the total 
(comoving) emissivity of the sources, ef , by 



N s 



- f 

Out 



± { JLy 
h P y v L ' 

3 + omg 
3a s 



. dv 



(1 



hpas 

T /iphs -1 Mpc~ 



(3) 



where As = 1.3 



+ (1.5,4.6) 
-(0.7,1.0) 



xl0 52 ,7 = -1.6-0.61og(As/1.3x 



10° 



as is the spectral index of the sources at frequencies 



above the photoelectric threshold (L v tx v~ as ), and the lcr 
and 2(j error ranges are provided. Here the approximation 
ef — eL has been made at the Lyman edge. This neglects 
the additional diffuse emission arising from the IGM itself 
which may enhance the total metagalactic photoionization 
rate over the direct rate from the sources by as much as an 
additional 50% (Meiksin & Madau 1993; Haardt & Madau 
1996). The source emissivity may thus be overestimated by 
this amount. In the spirit of providing an upper limit to the 
reionization redshift consistent with the measurements of 
the mean transmitted Lya flux, the more generous estimate 
of the source emissivity is adopted here. 

The evolution of asNs (comoving) is shown in Fig. 2 
(solid curve), using the fit to eL above. Allowing for a boost 
of 50% in the total emissivity by diffuse emission from the 
IGM reduces the required photon emission rate from sources 
by the same amount (long dashed line). These rates are 
compared with the minimum rate njrtjj 1 at which ionized 
photons must be produced per unit volume over a Hubble 
time, i H « (2/3)(H n 1 ^ 2 )- 1 (l + z)~ 3 ' 2 (for z » 1), to 
reionize the IGM. The two rates cross at z « 6 — 13, de- 
pending on the value of the spectral index. This coincidence 
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Figure 2. The evolution of the comoving production rate of ioniz- 
ing photons. The required rate inferred from measurements of the 
mean transmitted Lya flux is shown by the solid line. The long- 
dashed line shows the required source rate if diffuse recombination 
radiation from the IGM boosts the total metagalactic photoion- 
ization rate by 50%. Also shown are the minimal ionizing photon 
production rate njjtjj 1 required to ionize the IGM over a Hubble 
time (dotted line), and the rate of recombinations for clumping 
factors of C = 1 (dot-dashed line) and C = 5 (short dashed line). 
These latter three rates are shown for ag = 1. For alternative val- 
ues, they should be multiplied by the factor as- Also shown are 
the predicted contributions from QSO sources, assuming either a 
pure luminosity evolution model (filled squares) or pure density 
evolution model (open squares) for the QSO luminosity function. 
(See § 3.) 



in values is expected immediately following rcionization, and 
provides independent support (but not conclusive evidence) 
that reionization happened during this time. The smooth 
evolution of the required emissivity for z < 6 (Fig. 1), re- 
flected by the smoothness of the fit Ns curve in Fig. 2, sug- 
gests extrapolation to higher redshifts is reasonable. Doing 
so results in too low an emission rate of ionizing photons 
compared with the required minimum number at redshifts 
above z « 11 for as = 0.5, and above z m 8.5 for as > 1. 

For 2 < 15, the full recombination time at the average 
density of the IGM exceeds the Hubble time, and the effect 
of recombinations will be small. Because the IGM is non- 
uniform, however, the reionization rate will be much higher 
in clumpy regions, while smaller in underdense regions. The 
volume-averaged recombination time will then be reduced by 
the "clumping factor" C = (n H )/(nH) (Shapiro & Giroux 
1987; Meiksin & Madau 1993). Clumping will slow down the 
growth of ionization regions due to recombinations within 
them, requiring the recombined gas to be reionized again and 
so drawing on the overall ionizing photon budget. The value 
of the clumping factor is unknown. Numerical simulations 
prior to reionization suggest values of C ~ 30 at z = 8 — 10 
and C > 100 by z = 6 (Springel & Hernquist 2003). It is 



unclear, however, that these are the appropriate values to 
use. Once the gas is ionized, simulations suggest a clumping 
factor of C w 3 — 4 (at z — 3) (Sokasian, Abel & Hernquist 
2003). Of particular relevance is the number of photons per 
atom required to photoionize a clump (or minihalo) of gas, 
which requires a detailed numerical computation beyond the 
scope of current full cosmological simulations. The compu- 
tations of Shapiro, Iliev & Raga (2004) and Iliev, Shapiro 
& Raga (2004) have begun to examine this problem. They 
find an average of 2 — 5 ionizing photons are required per 
hydrogen atom, depending on redshift, with fewer photons 
required the harder the spectrum and the fainter the mean 
flux incident on a minihalo. 

In Fig. 2, the required minimum photon production 
rate, Cnut^,]., assuming the IGM is recombination domi- 
nated are shown with clumping factors of C = 1 and 5. 
Here, t Tec is the full hydrogen recombination time at the av- 
erage density of the IGM, t lac — l/n e a_g(T), where n e is 
the average electron number density in the IGM and as (T) 
is the radiative recombination rate to the excited levels of 
H i . (This will be taken at the reference temperature of 
T — 2 x 10 4 K for numerical values below.) For C = 5, 
the recombination time is shorter than the Hubble time for 
z > 4.4, and nearly half the Hubble time at z = 7. For 
this case, recombinations will play an important role in the 
reionization of the IGM. On the other hand, for C = 3 the 
recombination time is shorter than the Hubble time only for 
z > 6.6. The role recombinations played in the reionization 
of hydrogen in the IGM is therefore unclear. 

Madau, Haardt & Rees (1999) considered an IGM that 
was recombination dominated, with clumping factors of up 
to C — 30. Yan & Windhorst based their assessment for 
whether a class of objects could reionize the IGM on this 
value. Such a large value, however, is inconsistent with the 
estimated values of Ns from measurements of the mean 
transmitted Lya flux in Fig. 2. Scaling the curve for nnt^ 
(C — 1) up by a factor of 30 boosts the volume- averaged 
recombination rate to a value that exceeds the estimated 
upper limit on the source rate by a factor of 3 at z = 6, 
even allowing for a hard source spectrum of as = 0.5. (The 
discrepancy is even larger for softer sources.) This suggests 
C < 10 is a more realistic upper limit. (The limit is even 
more restrictive for softer sources.) 

The epoch of reionization is quantified using the ion- 
ization porosity parameter. The H n ionization porosity, 
Qhii, evolves according to (Meiksin & Madau 1993; Madau 
et al. 1999): 



dz 



N s {z) C{z)Qi 



n H (0) t lec (z) 



dt_ 

dz ' 



(4) 



where n H (0) is the comoving number density of hydrogen 
atoms. If the recombination time is long compared with the 
Hubble time, the recombination term in Eq. (4) may be 
neglected. Without a detailed knowledge of the clumping 
factor and its evolution, including feedback effects, a com- 
putation still beyond the capability of current large-scale 
gravity+hydrodynamics+radiative transfer simulations, the 
role of recombinations cannot be determined in detail. Here, 
C = 1 will generally be assumed, so that only an upper red- 
shift limit to the epoch of reionization is estimated. Signif- 
icant clumping will serve to further delay complete reion- 
ization (Q > 1). The effect of clumping on slowing down 
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the reionization process may be particularly important for 
He n reionization (Madau & Meiksin 1994). 

Using the expression Eq. (3) in Eq. (4) and integrating 
gives, in the absence of recombinations, 



Qhii(z) = 



N s (z)t H (z) 
n H (0)(l-§ 7 ) 



/ 1 + z 
U + 



(5) 



where Zi is the turn-on redshift of the sources. In the 
presence of recombinations (and for a redshift independent 
clumping factor and IGM temperature), the solution is 



Qmi(z) = 



Ns(z)tn(z) i- 
n H (0) Th 

/■th(z») 

x / dre~ 

Jtk(z) 



\z)e 



(6) 



where th(z) is the Hubble time measured in units of the re- 
combination time (including the clumping factor), Tn(z) = 
Ctn(z) /t rcc (z) ~ 0.016C(1 + z) 3/2 , and is equivalent to the 
mean number of recombinations in the IGM over a Hubble 
time. Provided 2 — 27/3 is not a positive integer, the inte- 
gral may be expressed in terms of the incomplete gamma 
function F(a,x) = J*J° dt e - t t a - 1 . For 7 > -3 (7 ^ 0, ±§), 
and Zi — > 00, 



Qmi(z) = 



Ns(z)M z) 
n H (0)(l 



§7) 



1 4. 3th 1 3r H 

2 7 2 7 (1 + | 7 ) 



m - e 



(2 + §7,1*) 



(7) 



When 2 — 27/3 is a postive integer, the solution is 

VHIl(z) - 7-T e 'iV2-2 7 /3(TH), (8) 

nH(Uj 

where E n {x) — dte~ xt t~ n is an exponential integral. For 
the special case of a constant comoving source emissivity 
(7 = 0), Eq. 8 recovers eq. (22) of Madau et al. (1999). 
In the absence of recombinations, th — » and Eq. 5 (for 
Zi >> z) is recovered. The asymptotic series representation 
of Eq. 6 for th(z) > 1 (and Zi >> z) is 



Qhii(z) ~ C- 1 



N S (z)t ICC {z) 
n H (0) 



S7j 



r£(*) ' 



(9) 



where (o)„ = T(a + n)/T(a) = a(a + l)...(a + n - 1). The 
first term (n = 0) is the approximation eq. (23) of Madau 
ct al. (1999), valid in the limit th >> 1. Similar expressions 
follow for the He in porosity, QHein, if h-h(O) is replaced 
by riHe(O), and Ns and £ rec refer to the production rate of 
He ii ionizing photons and the recombination rate to He n . 

The key assumption is now made that the sources of 
ionization at z < 6 are drawn from the same population 
as that responsible for the reionization of the IGM. Under 
this assumption, the porosity Qhii may be estimated from 
Eq. (3) for Ns, extrapolated to z > 6. This seems a rea- 
sonable assumption to make as the inferred emissivity in 
Fig. 1 evolves smoothly over 3 < z < 6 in spite of the strong 
downward turn in the measured mean Lya transmitted flux 
at z > 5. To make the estimate, however, it is first necessary 
to be more specific concerning the factor (3 + omgJ/qs in 
Eq. (3). 

When QSOs are the sources of reionization, the UV 
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Figure 3. The evolution of the effective spectral index a e g (see 
text) for a starburst of solar mctallicity with a constant star for- 
mation rate (solid line) or a sudden burst of star formation (short 
dashed line). Also shown is the evolution assuming a mctallicity 
Z = 0.05Zq and a constant star formation rate (long dashed 
line). 



metagalactic background hardens at high redshifts beyond 
the input spectral shape because high energy photons travel 
further than photons near the photoelectric threshold before 
being absorbed by the IGM. Here o?mg ~ 1 is adopted based 
on the results of Haardt & Madau, including re-emission 
from the IGM. Although the spectral index should be com- 
puted self-consistently from the models considered here, the 
uncertainty in this value introduces only a small uncertainty 
in the total photoionization rate. It is noted that o?mg may 
be double this if young star forming galaxies with a Salpeter 
IMF and solar abundances are the sources of ionization, as 
they have much softer spectra than QSOs. 

More critical is the value of as, which sets the to- 
tal number of ionizing photons released by a source in 
terms of its specific luminosity at the photoelectric thresh- 
old. The value for a QSO is uncertain. Results based on 
Hubble Space Telescope (HST) data give as = 1.76 ± 0.12 
(Telfer et al. 2002). More recent results based on Far Ul- 
traviolet Spectroscopic Explorer (FUSE) data give instead 
a s = 0.56±g;§| (Scott et al. 2004). The origin of the dis- 
crepancy is unclear. The FUSE sample is at lower redshift 
and the sources have lower luminosities, so the difference 
may reflect a redshift and/or luminosity dependence. The 
corrections due to intervening absorption by the IGM, how- 
ever, are significantly smaller for the FUSE sample than for 
the HST sample, so that the FUSE result may be a more 
reliable estimate of the spectral index of the mean QSO 
spectrum. 

For a non-power-law source like a galaxy, an effective 
spectral index a| ff may be defined according to 
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„ cff — , in , 

as = — 7~, (10) 
hpN s 

where Ll is the specific luminosity of the source at the Ly- 
man edge. The effective spectral index is shown in Fig. 3 for 
starburst galaxies for both a fading sudden burst and for a 
constant star formation rate, using Starburst99 (Leitherer 
et al. 1999), updated with the stellar spectra of young stars 
from Smith, Norris & Crowther (2002). A Salpeter IMF is 
assumed and solar abundances. While the spectral index of 
a fading burst quickly climbs to values exceeding 10 (after 
4 x 10 7 yrs), in the case of continuous star formation, the ef- 
fective spectral index instead reaches an asymptotic value of 
Qj 8 « 2.3. For continuous star formation with Population 
II abundances (Z = 0.05Zq), the asymptotic value of the 
effective spectral index is a| ff 1.8. This is comparable to 
the effective spectral index for a blackbody spectrum with 
a temperature of T — 5 x 10 4 K , used by Stiavelli et al. to 
represent a galaxy dominated by Population II stars. They 
represented Population III stars (essentially no metals) by a 
blackbody with T = 10 5 K , for which a cB « 0.5. 

One implication of the rapidly climbing value of the ef- 
fective index in the case of a fading burst is that the ionizing 
efficiency of the galaxy declines rapidly. If star formation in 
a bursting galaxy were episodic, with a duty cycle of 50% 
and a cycle time longer than 2 x 10 7 yrs, the net contribu- 
tion of the population to the number of ionizing photons is 
reduced by 50% from the case of continuous star formation. 
This complication is not considered here. 

Because of the uncertainty in the spectral index, three 
values are considered here for the spectrum near the hy- 
drogen Lyman edge: a hard spectrum with as = 0.5 corre- 
sponding to a hard spectrum QSO or Population III dom- 
inated galaxy, as = 1.8 corresponding to a soft spectrum 
QSO or Population II dominated galaxy, and as = 2.3 cor- 
responding to a starburst with solar metallicity. 

In Fig. 4, the evolution of the porosity parameter is 
shown for hard and soft source spectra, and presuming the 
sources turned on at very high redshifts (zi — > oo). Within 
the observational uncertainties, the earliest epoch of reion- 
ization for hard sources (as = 0.5), like QSOs or Pop III 
stars, is z r i < 11 (2a upper limit). For soft sources (a = 2.3) 
like starburst galaxies, z ri < 5.5 (2a). Intermediate sources 
with as = 1.8, like soft spectrum QSOs or Pop II stars, may 
ionize the IGM only at z r i < 6 (2a). While hard spectrum 
sources are well able to photoionize the IGM prior to z = 6, 
reionization by softer spectrum sources like soft spectrum 
QSOs or starbursts (unless dominated by Pop III stars), is 
more problematic. 

The estimate provided neglects several factors which 
will further delay the epoch of reionization. In particular, 
reionization will be delayed if the galaxies turn on later (fi- 
nite Zi) or if the clumping factor is large. The source emissiv- 
ity at the Lyman edge was also equated to the total required 
metagalactic emissivity, neglecting the additional ionizing 
photons arising from the IGM itself. Taking this into account 
could reduce the estimated source emissivity by as much as 
50% (see above), and so lower the overall predicted rate of 
ionizing photons by this same amount, delaying reionization 
even further. 

The 2a upper limits on z r i for soft and intermediate 
spectrum sources lie well below the WMAP la and 2a lower 
limits of z r i > 14 and 11, respectively (Kogut et al. ). The 
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Figure 4. The evolution of the ionization porosity parameter 
Qhii f° r hydrogen reionization based on the metagalactic ionizing 
emissivity estimated from measurements of the mean transmitted 
Lyct flux. The solid lines show the expected evolution assuming 
a source spectral index of, from right to left, as = 0.5, 1.8 and 
2.3. The dashed lines show the corresponding evolution for the 
la upper limits on the estimated emissivity. A clumping factor 
of C = 1 is assumed for the above curves. Also shown are the 
corresponding 2cr upper limits for C = 3 (dot-dashed curves). 



prediction for as = 0.5 is just barely compatible with the 
WMAP 3a lower limit of z r i > 8. It should be noted that 
Spergel et al. (2003) find the somewhat broader reionization 
redshift constraint z T \ = 17 ± 4 after combining with other 
data sets. This still disagrees with the intermediate and soft 
spectra predictions at the 3<r level, but agrees with the pre- 
diction for as = 0.5 at the Si 2a level. The results suggest 
that if the ionization of the IGM at z < 6 is provided by 
soft or intermediate spectrum sources, like starbursts dom- 
inated by Pop II stars, they are not members of the same 
population that reionized the IGM, unless the population as 
a whole was a much more prodigious source of ionizing radi- 
ation at earlier times (z > 6). Sources with a hard spectrum 
(as = 0.5), like Pop III stars or hard spectrum QSOs, are 
preferred. 

Gnedin (2004) similarly found that for his numerical 
simulations to match the mean Lya transmitted flux mea- 
surements at z < 6, the epoch of reionization (defined some- 
what differently in terms of the volume- weighted mean IGM 
ionization fraction), must have occurred at z r i ~ 6. The 
results here are much more generic and not dependent on 
any particular prescription for cosmic star formation, as in 
Gnedin's simulations. Gnedin finds that the predicted reion- 
ization epoch is determined by a single parameter, the ultra- 
violet emissivity parameter. Setting this parameter to match 
the mean Lya optical depth measurements at z < 6 requires 
reionization to have occurred at z « 6. The one parameter 
nature of Gnedin's models appears to be a consequence of 
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the assumed nature of the ionizing sources formed as part 
of the growth of cosmological structure. The analysis here 
shows that reionization is determined by three parameters, 
expressed by As, 7 and as. The constraint on the reioniza- 
tion epoch found here follows only from the assumption that 
the population of sources which dominates the ionization 
background at z < 6 is the same population that reionized 
the IGM. The Lya optical depth measurements at z < 6 
are used to determine As and 7, leaving as free to set the 
reionization epoch. 

In the next section, the QSO contributions to the meta- 
galactic ionizing background and to the total budget of ion- 
izing photons are estimated based on the QSO counts in 
recent surveys. 



3 THE CONTRIBUTION OF QSOS TO THE 
IONIZING PHOTON BUDGET 

The QSO luminosity function for z < 2 is well- 
established. The QSO population detected by the 2dF 
(Croom et al. 2004) is well-fit at least up to z < 2.1 by 
a pure luminosity evolution (PLE) model with the double 
power- law form (Boyle 1988): 

4>( M ,z) = — 4(1 „ /3l)[M _ M » (z)] - 10 0.4(1-/3 2 )[M-M*( Z )] '( n ) 

where <f>* is a fixed spatial normalization factor and M*(z) 
is an evolving break magnitude. 

For z > 3, the model parameters are less clear since 
surveys have generally only detected QSOs brighter than 
the break magnitude. The most complete of these surveys 
is the QSO survey at z > 3.6 carried out as part of the 
SDSS (Fan et al. 2001, 2004). An exception is the survey of 
Steidel et al. (2002), in which Active Galactic Nuclei (AGN) 
dimmer than the break magnitude were discovered as part 
of a Lyman break galaxy survey at z ~ 3. 

In this section, the contribution of QSOs to the ioniza- 
tion rate of the IGM over the redshift range 3 < z < 6 is 
estimated based on the findings of Fan et al. (2001, 2004), 
and as also constrained by the low luminosity AGN data of 
Steidel et al. (2002). This is done for both pure luminos- 
ity evolution and pure density evolution (PDE) models, for 
which M* is fixed in Eq. (11) but cj>* is allowed to evolve, as 
an attempt to bracket the uncertainty in the evolution. 

The luminosity function parameters in Eq. (11) are es- 
timated using a maximum likelihood procedure. For any 
given set of parameters, the predicted numbers of QSOs in 
a set of absolute magnitude and redshift bins are given by 
(J-ij — <P{Mi, Zj)V^ tl (M i )dM i for magnitude bin i and redshift 
bin j, where V^ s (Mi) is the effective survey volume probed 
for QSOs with absolute magnitudes within the bin centred 
on Mi and redshifts in the bin centred on Zj. All absolute 
magnitudes are taken at the published values at 1450A, and 
a flat cosmological model with Qm = 0.35 and h = 0.65 
is used for the fits for consistency with Fan et al. (2001, 
2004). (The magnitudes and effective volumes from Steidel 
et al. (2002) and Hunt et al. are transformed to this cosmol- 
ogy.) The likelihood of the model compared with the actual 
numbers Nij of QSOs detected is then given by 




M-M* 

Figure 5. The maximum likelihood fits to the QSO luminosity 
function for the PLE (solid line) and PDE (dashed line) models. 
The measured luminosity function values have been renormalizcd 
by <f>* and M* for the PLE model (solid points) and the PDE 
model (open points). The error bars are Poisson, based on the 
maximum likelihood fits. 



where Poisson probabilities are required because of the low 
number of objects per bin. 

The likelihoods are computed on a multi-dimensional 
(linear) grid in /3i, fii, 4>* and M* . (This corresponds to 
assuming a uniform prior for each parameter.) Here f}\ is 
assumed to correspond to the low luminosity end of the QSO 
luminosity function and /?2 to the high end. For the PLE 
models, the same <f>* is used for all redshift bins, but M* is 
stepped over separately for each redshift bin. Conversely, for 
the PDE models the same M* is adopted for all redshifts, 
but <j)* is allowed to vary for each redshift bin. Over 10 9 
models are constructed for the PLE case and 10 s for the 
PDE case. 

The parameter values for the maximum likelihood fit of 
the PLE model are (3i = 1.24, /3 2 = 2.70 and <j>* = 3.25 x 
10~ 7 Mpc -3 mag -1 (comoving). The values for M*(z) are 
given in Table 1. The parameters for the PDE fit are /3i — 
1.12, #2 = 2.90 and M* = -25.60. The values for </>*(z) are 
given in Table 1. 

The goodness-of-fit of the maximum likelihood models 
is estimated using the x 2 -test for the predicted vs detected 
number of QSO sources for all 15 available bins in redshift 
and magnitude (one of which is a cumulative magnitude 
bin), and allowing for 7 degrees of freedom (after accounting 
for the 8 model parameters). For the PLE model, \ 2 = 6.8 
and p(x 2 > 6.8) = 0.45. For the PDE model, \ 2 = 8.5 and 
p(x 2 > 8.5) = 0.29. Both fits are acceptable. The maximum 
likelihood fits are shown in Fig. 5, along with the (renormal- 
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ized) measured values of the luminosity function and their 
Poisson errors, which are based on the fit models. 

The maximum likelihood parameter values, each 
marginalised over all other model parameters, and with 
their equivalent la error ranges, are, for the PLE model, 
Pi_ = 1.24±g;i5, (3 2 = 2.2±° 3 2 and 0* = 2.8^;° 8 x 
10 7 Mpc 3 mag 1 (comoving). The marginalised values for 
M*(z) are given in Table 1. For the PDE model, the val- 
ues are j3i < 1.1, (3 2 = 2-Sto.t and M * = — 25.4±g;|. The 
marginalised values for <j>*(z) are given in Table 1. The upper 
limit on (5\ results since smaller values were not considered. 

The results on f5\ and 02 agree well with the findings of 
Hunt et al. and Fan et al. (2001, 2004). Hunt et al. obtain 
0i = 1.24 ± 0.07 at z « 3. Fan et al. (2001) obtain, for 
3.6 < 2 < 5.0, (3 2 = 2.58 ± 0.23, while at z > 5.7, Fan 
et al. (2004) obtain f3 2 = 3.2±° *. (All errors are Iff.) The 
PLE model estimate of f} 2 found here disagrees with the 
z > 5.7 value at the 98% confidence level, suggesting some 
departure from PLE, but the evidence is not conclusive. 

For 0.4 < z < 2.1, Croom et al. obtain /3i = 1.0 ± 0.1 
and /?2 = 3.25 ± 0.05 for a PLE model with a break mag- 
nitude varying as an exponential of lookback time. (The re- 
sults for a quadratic redshift dependence are similar.) There 
is again only marginal agreement with /3 2 . There is greater 
disagreement over the normalizations. Croom et al. report 
<)>* B = 1.8 x 10~ 6 Mpc~ 3 mag -1 for the B-band luminosity 
function. This value (including an adjustment for the small 
differences in assumed cosmologies), exceeds the values for 
<j>* found for the PLE model here by a factor of 6lg , signif- 
icant at the 5a level. It should, however, be noted that the 
normalizations may only be fairly compared provided the 
ratio of B— band to 1450A luminosities was itself luminosity- 
independent. A steep correlation between spectral slope and 
B— band magnitude for QSOs could in principle account for 
the difference in normalizations. Currently there is no ev- 
idence for a correlation sufficiently steep. It thus appears 
there has been substantial density evolution from z ~ 3 to 
z<2. 

For each set of luminosity function parameters, the inte- 
grated emissivity values at the Lyman edge for the PLE and 
PDE models are computed from e L = J <j>(M)L L (M)dM , 
where Ll{M) is the frequency specific QSO luminosity at 
the Lyman edge corresponding to a given absolute magni- 
tude M at 1450A, and assuming L v tx v~ QS ' between 1450A 
and the Lyman break. The resulting marginalised maximum 
likelihood emissivity values (transformed to a flat cosmol- 
ogy with Q.m = 0.3 and h — 0.7), are shown in Fig. 1 
for the PLE model (filled squares) and the PDE model 
(open squares), along with their la errors. The correspond- 
ing marginalised estimates for the ionizing photon produc- 
tion rates [asNs = th/hp) are shown in Fig. 2. 

As shown in Fig. 1, for z £ 3, all of the required emis- 
sivity is accounted for by QSOs. At z > 4.5, however, it 
appears QSOs are inadequate sources of ionizing photons, 
falling short by factors of at least 2 — 3, depending on the 
amount of diffuse radiation. Additional sources are required 
to make up the difference. Within the uncertainties, the QSO 
deficit could be as high as a factor of 30 by z « 5.5. 

The deficit of ionizing photons from QSOs is apparent 
in Fig. 2 as well. For z > 5.5, QSOs fall short of the ion- 
izing photon production rate required for consistency with 
the measured mean transmitted Lya flux values (solid and 
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Figure 6. The probability distributions for the epoch of rcioniza- 
tion for hydrogen (solid lines) and helium (dashed lines) assuming 
a QSO spectral index of as = 1.8 (heavy lines) or 0.5 (light lines). 

long-dashed lines) by at least a factor of 5. The la upper 
limit on the QSO rate for the PLE case at z = 5.7, how- 
ever, comes close to the bare minimum rate required to ion- 
ize all the baryons over a Hubble time, assuming a hard 
(as — 0.5) source spectrum. While QSOs could not have 
ionized the IGM acting alone, they still may have provided 
a non-negligible contribution of ionizing photons during the 
reionization epoch. 



4 DISCUSSION 

Any assessment of the contribution of a population of 
sources to the ionizing photon budget of the IGM depends 
on several unknown factors. Principal among these are the 
spectral index ckmg of the UV background, the contribution 
of Lyman Limit Systems to the attenuation length at the Ly- 
man edge, the amount of diffuse emission arising from the 
IGM, and the clumping factor of the IGM. In this section, 
the role each of these factors may play is discussed. 



4.1 The sources of ionization of the IGM 

Comparison of the results of numerical simulations of the 
IGM with measurements of the mean transmitted Lya flux 
through the IGM constrains the ionization rate per hydro- 
gen atom (or per He n atom for He n Lya ). This is related 
to the emissivity by a factor of 3 + omg (Eq. 1). The re- 
sults of Haardt & Madau suggest qmg ~ 1 shortward of the 
hydrogen Lyman edge, but this is due in large part to the 
contribution of redshifted He n recombinant Lya and two- 
photon emission. If most of the helium is still in the form 
of He ii , the spectrum may be harder, possibly even with 
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omg < 0. In this case, the required emissivity could be re- 
duced by as much as a factor of two, bringing it much more 
in line with the expected QSO emissivity. 

Meiksin & White (2004) appealed to the uncertainty 
in the values of omg and the QSO break magnitude (M*), 
adjusting these to allow QSOs to dominate the UV ioniz- 
ing background at high redshifts. The assessment here does 
not rule out this possibility. Since QSOs acting alone would 
not be able to photoionize the IGM by z = 6, however, it 
seems likely the UV background will have a large contri- 
bution from other sources, such as galaxies. (Although this 
need not be the case if reionization occurred at z » 6.) 
While this would suppress the correlations in the UV back- 
ground predicted by Meiksin & White (2004), and the pre- 
dicted enhancement of the Lya flux correlations in the IGM, 
the suppression may be by no more than a factor of a few. In 
this case, the contribution of QSOs to the UV background 
may still be detectable through the fluctuations they induce. 

The ionization rate predicted for a given emissivity also 
depends on the attenuation length ro- The values used here 
were computed self-consistently from numerical simulations, 
but these simulations do not include radiative transfer and 
therefore are not able to accurately predict the contribu- 
tion of Lyman Limit Systems. These are accounted for by 
direct observations for z £ 4, but their abundances are not 
known at higher redshifts. The attenuation length may be 
~ 30% smaller than the values adopted here, resulting in 
an increase in the required emissivity by the same amount. 
This would result in a further reduction in the fraction of 
the ionizing background provided by QSOs. A similar effect 
may arise if the emissivity of the sources rapidly declines 
with redshift, in which case the effective path length over 
which the emissivity is integrated to obtain the ionization 
rate may be comparable or even shorter than the attenua- 
tion length, which also would permit a higher emissivity at 
a given redshift (Meiksin & White 2003). 

If a substantial fraction of the required emissivity 
(~ 50%) arises from diffuse recombination radiation within 
the IGM, then the requirements for the source emissivity are 
reduced by the same amount. The fractional contribution of 
QSOs to the required source emissivity would then increase. 
The difference between the required source emissivity and 
the predicted emissivity from QSOs for the PLE and PDE 
cases is shown in Fig. 1, assuming the diffuse emission from 
the IGM contributes at 50% the source level. QSOs pro- 
vide more than half of the required source emissivity up to 
z 4.5, and less beyond. The curves showing the difference 
describe the emissivity contributed by other sources, pre- 
sumably galaxies. 

4.2 The sources of reionization of the IGM 

The reionization redshifts estimated here rest on the as- 
sumption that the population of sources which dominates 
the UV ionizing background at z < 6 is the same popu- 
lation responsible for reionization. This is to be expected 
if reionization occurred near z « 6, as argued by Becker 
et al. , but there is nothing to preclude a much earlier reion- 
ization epoch (z >> 6), or even a series of several reioniza- 
tion epochs followed by recombination. On the other hand, 
the mild evolution in the emissivity required to match the 
measured Lya optical depths at z < 6 suggests there was 



no abrupt change in the nature of the sources dominating 
the UV background at z ~ 6, so that the emissivity may be 
reasonably extrapolated to higher redshifts. 

The estimated reionization redshifts were all upper lim- 
its, neglecting any excess clumping (C > 1) in the reion- 
ized gas. If the clumping factor is moderate (C ft* 2 — 4), as 
suggested by recent numerical computations, then recombi- 
nations may not delay the reionization epoch very signifi- 
cantly. The evolution of the porosity for C — 3, assuming 
the 2cr upper limit on the emissivity, is shown in Fig. 4. For 
as — 2.3 and 1.8, reionization now occurs well below z = 6. 
For as = 0.5, the 2a upper limit of z r i < 9 is still compat- 
ible with the SDSS data (Becker et al. ), but barely with 
the WMAP 2(j lower limit of Spergel et al. The expected 
value of the reionization redshift for as = 0.5 and C — 3 is 
z T i = 7.1, which is compatible with the WMAP 2.5a lower 
limit of Spergel et al. For higher values of the clumping fac- 
tor (C ~ 5 — 10), reionization for the various spectral index 
cases considered here would occur substantially later. 

If the reionization of minihalos doubles the reionization 
requirements (Shapiro et al. 2004, Iliev et al. ), then reion- 
ization for as = 0.5 is not expected to occur until z r i ^ 6.2. 
While this is consistent with the data of Becker et al. , it 
is consistent with the WMAP lower limit on z r i of Spergel 
et al. at only the 2.7a level. 

The estimates here assumed the sources turned on at in- 
finitely high redshifts (zi >> z T i). Reionization would occur 
somewhat later for finite values of Zi. 

The estimates for the reionization epoch were also based 
on the assumption that all the required emissivity was pro- 
vided by the sources, while in fact diffuse emission from the 
IGM could contribute an additional ~ 50%. In this case, the 
expected source contribution would be reduced by 50%, and 
(1 + Zri) reduced by about 15%. For as = 0.5, this would 
push the 2a upper limit on z r \ down from 11 to 9. 

Reionization may have occurred earlier if the emissivity 
values used here are underestimated, as would be the case 
if the attenuation length of Lyman photons were shorter 
than the values used here (see § 4.1 above). Even a doubling 
of the emissivities, however, would increase the expected 
values of (1 + z r \) by only about 20%, or z-a ~ 10 for as = 
0.5. For as = 1.8 and 2.3, the expected redshifts for the 
reionization epoch may be increased to z T \ « 6, still too 
low to match the WMAP constraints, but consistent with 
the SDSS Lya optical depth measurements. At much higher 
redshifts, the reionization epoch will be limited by the effect 
of recombinations even for small clumping factors. 

In the scenarios above, the expected redshift of 
He ii reionization is always above z r i > 11 for as = 0.5 
(even for C = 3), so that helium may have been ionized 
very early if the source spectra continue as hard power laws 
shortward of the He n Lyman edge. If the spectrum softens 
to as > 1.8 shortward of the He n Lyman edge, however, 
He ii reionization would not have occurred prior to z r i £ 6 
(except if the total emissivity is doubled), but generally prior 
to 2 r i ^ 3, unless He n ionizing photons are produced at a 
negligible rate, as is expected if the ionizing sources are star 
forming galaxies (Smith et al. ). 

QSOs are found to provide too few ionizing photons 
at z > 5 to reionize the hydrogen in the IGM. This may 
be quantified in terms of the ionization porosity. Because 
of the rapid decline in the photon production rate, Eq. (6) 
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no longer applies. Instead, the decline in Fig. 2 is better 
modelled as an exponential in redshift: asNs = AseT' ,z ■ 
The solution of Eq. (4) (with recombinations neglected) is 

Qmi= n H (o) [ l ~ 2 y +2e y r {2' y )\' (13) 

where y = 7(1 + z), and z t — ► 00 is assumed. A similar 
expression follows for QhcIii on replacing nn(0) by riHe(O) 
and N s by iV s /4 as . 

Because of the broad, non-gaussian probability distribu- 
tion for Ns from QSOs, the probability distribution for the 
reionization redshift is determined from a maximum likeli- 
hood estimate using the probability distributions found in 
§ 3 for the QSO emissivity. For the PLE case, the parameters 
of the maximum likelihood model are As = 5.7 x 10 5 ph s~ 
and 7 = 1.26. For the PDE case, A s = 6.1 x lCPphs -1 
and 7 = 1.37. The marginalised maximum likelihood values 
for the parameters are, for the PLE case, As = 5.2I3' 2 x 
10 52 phs _1 and 7 = 1.65 + 0.10 log(As/5.2 x 10 62 ). For 
the PDE case, A s = 5.7± 7 3 t x lO^phs" 1 and 7 = 1.74 + 
0.121og(4s/5.7 x 10 82 ). 

The resulting probability distributions for the reioniza- 
tion redshifts, both for H 11 and He in , are shown for the 
PLE and PDE cases in Fig. 6. Results for both as = 1.8 and 
0.5 are provided. For as = 1.8, the hydrogen and helium 
reionization redshifts nearly coincide (since the coefficients 
for Qhii and QhcIii are nearly identical for this value of 
the spectral index) , and correspond to reionization redshifts 
z r i ~ 3. While inconsistent with the H 1 data, measurements 
of the mean transmitted He 11 Lya flux reveal the absence 
of a measurable flux by z ~ 3.5 (Zheng et al. 2003). Any fu- 
ture detections of a finite He n Lya optical depth at z > 3.5 
would suggest the QSO spectra are harder than as = 1.8 
shortward of the He n Lyman edge. For as = 0.5, helium 
reionization may occur as early as 2 fx 5 by QSOs alone. 
By contrast, even for this hard a spectrum, QSOs alone are 
unable to reionize the hydrogen in the IGM any earlier than 
z * 4. 

If helium was not fully ionized until z « 3, then the 
UV background may have undergone a second change in its 
spectral shape shortward of the hydrogen Lyman edge due 
to the release of (redshifted) recombinant He n photons and 
two-photon emission (Haardt & Madau). (The first change 
may have been at z < 4.5 when QSOs began to domi- 
nate the source emissivity.) There would of course also be a 
marked change shortward of the He 11 Lyman edge. In ad- 
dition to the indirect effect on the temperature of the IGM, 
these changes may be detectable through their effects on the 
relative abundances of ionization states of metal ions (eg, 
Si iv /C iv ), for which there is some evidence at z « 3.1 
(Songaila & Cowie 1996). 



4.3 Did low luminosity AGN reionize the IGM? 

The implication that the reionization sources had hard spec- 
tra (as — 0.5) would suggest that star forming galaxies 
(or star clusters) dominated by Pop HI stars are the only 
viable candidates for reionizing the IGM. QSOs, the only 
other hard spectra sources known, decline too rapidly in 
their number at high redshifts. Reionization soley by Pop III 
stars, however, has its own difficulties, as the stars would 



rapidly pollute their environment with metals, so that reion- 
ization would be accomplished predominantly by Pop II 
stars (Schneider et al. 2002; Ricotti & Ostriker 2004). As 
shown above, Pop II stars are just barely able to reionize the 
IGM by z — 6 without violating the emissivity requirements 
needed to match the measurements of the mean transmitted 
Lya flux at z < 6. They fall well short of the WMAP reion- 
ization requirement. Pop III stars (or Pop II stars) could 
only have reionized the IGM if their emissivity rapidly de- 
clined from the time of reionization to z < 6. 

Madau et al. (2004) discuss miniquasars as an alter- 
native reionization source, arguing that they may arise nat- 
urally in the same systems forming Pop III stars through 
accretion onto black holes formed from collapsed massive 
stars. Dijkstra, Haiman & Loeb (2004), however, claim that 
this and any scenario in which reionization was performed 
by quasar-like objects are excluded by constraints from the 
soft X-ray background. In this section, the possibility that 
low luminosity AGN reionized the IGM is re-explored. 

From an energetic perspective, the requirements that 
galaxies harbouring active nuclei at high redshifts must 
satisfy for reionization are readily met. A black hole of 
mass 10 6 ' 5 Mq would shine with an Eddington luminosity of 
L E dd — 5 x 10 44 ergss -1 , producing a specific luminosity at 
the Lyman edge of Ll ~ 0.1L E dd/^L ~ 10 28 ergss -1 Hz~ x 
(allowing for a contribution near the Lyman edge to the 
bolometric luminosity of 10%). At z = 6, Eq. 2 gives a total 
(comoving) emissivity of el ~ 3 x 10 24 ergss - Hz -1 Mpc~ 3 . 
The corresponding comoving space density of galaxies har- 
bouring such black holes is then n 9 ~ 2x 10~ 4 Mpc~ 3 , only a 
small percentage of galaxies today. The magnitude at 1450A 
would then be M1450 ~ —19 (assuming /„ ~ v~ ). This 
corresponds to a z' magnitude of £ 27 at z ~ 6 and ^ 28 
at z ~ 8. The number of sources per unit redshift per square 
arcminute at 6 < z < 8 would then be dN/dz « 1/2. It is 
possible that as many as a dozen of such objects appear as 
i'-dropouts in the UDF (see, eg, Bunker et al. 2004). At such 
low luminosities, however, they would have escaped detec- 
tion in the GOODS fields and in the z « 3 survey of Steidel 
et al. (2002). It does imply a sharp turn up at faint mag- 
nitudes in the number of sources compared with the Hunt 
ct al. luminosity function at M1450 = —20. The comoving 
density of the objects, however, may have declined substan- 
tially by z w 3, as they are no longer required to maintain 
the ionization of the IGM, since the known QSOs are ade- 
quate (Fig. 1). The excess towards fainter magnitudes may 
then be only mild. 

The low luminosities and abundance of the objects also 
place a minimal demand on the number of ionizing pho- 
tons required per hydrogen atom to reionize the IGM. The 
typical (proper) distance between the sources at z = 11 is 
~ 1 Mpc, corresponding to a flux of ionizing photons of 
F ~ 10 53 phs" 1 Mpc -2 , or F = 0.01 - 0.1 in the units of 
Iliev et al. Such low fluxes require only 1-2 ionizing photons 
per hydrogen atom to photoionize a minihalo (Iliev et al. ) 
as a consequence of photoevaporation. 

The objects would also be consistent with constraints 
imposed by the soft X-ray background. The estimate of Di- 
jkstra et al. was based on assuming a mean requirement of 
10 ionizing photons per hydrogen atom to reionize the IGM. 
They then argued that any quasar-like sources which pro- 
duced sufficient photons to reionize the IGM would exceed 
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by a factor of about 4-6 the component of the residual soft 
X-ray background unaccounted for by discrete sources. As 
this factor scales linearly with the number of ionizing pho- 
tons required per atom, a requirement of only 2 photons per 
atom would bring the level into agreement with the residual 
background. (This argument may apply to miniquasars as 
well.) 

Even without appealing to the decrease in the number 
of required ionizing photons per atom, however, the soft X- 
ray background constraint may be avoided because of the 
presumed hard spectrum (as = 0.5) of the sources. Dijkstra 
ct al. based their estimates for QSOs on the spectrum of 
Sazonov, Ostriker & Sunyaev (2004), who assumed a mean 
QSO energy spectrum varying as (their eq. [14]) (F(E)) ~ 
E-°- 6 for leV < E < 10 eV, and (F(E)) ~ E' 1 -'' ' e B / 2keV 
for 10 cV < E < 2keV, based on the HST data of Telfer 
et al. Adopting as = 0.5 at E > 10 eV, however, can in- 
crease the number of ionizing photons while leaving the X- 
ray constraints unchanged. For instance, adopting the form 
(with E in units of keV) 

/ p* v _ n / 0.63£-°' 6 0.001 < E < 0.01, , 

preserves the optical to X-ray spectral index a ox between 
2500A and 2 keV (as well as the normalization at 2 keV), 
but produces 3 times as many hydrogen ionizing photons 
as the corresponding model of Sazonov et al. (It also pre- 
dicts Ll — 0.11/boi/^L, used above to relate Ll to Z/Edd-) 
This again would bring the predicted contribution of the 
sources to the soft X-ray background into agreement with 
the residual component. It thus appears that hard spectrum 
quasar-like objects may reionize the IGM without exceeding 
the limits of the soft X-ray background limits. 



5 CONCLUSIONS 

A relation between the mean transmitted Lya flux and ion- 
ization porosity of the IGM is exploited to obtain generic 
limits on the epoch of reionization. Under the assumption 
that the sources which dominate the UV ionization back- 
ground at z < 6 are members of the same population of 
sources that reionized the IGM, it is shown that the epoch of 
reionization must have occurred at z r i < 11 (2cr upper limit) 
for hard spectrum (as = 0.5) sources like QSOs, AGN and 
Pop III dominated star forming galaxies, z rl < 6 (2a) for 
moderate spectrum (as = 1.8) sources like soft spectrum 
QSOs, AGN and Pop II dominated star forming galaxies, 
and 2 r i < 5.5 (2a) for soft spectrum (as = 2.3) sources like 
star forming galaxies with solar abundances. 

After reionization, it is possible that no single popula- 
tion of sources dominated the UV ionizing background down 
to 2 w 2. Instead the UV ionizing background may have un- 
dergone multiple qualitative changes in its spectral shape 
as the nature of the dominant ionizing sources changed. 
An assessment of the contribution of QSOs to the budget 
of ionizing photons based on recent high redshift QSO sur- 
veys suggest that QSOs alone are adequate for providing the 
ionization rate required to reproduce the mean Lya optical 
depth measurements at z £ 3. QSOs produce at least half of 
the required emissivity at z « 3.5, but fall short by a factor 
of at least 2 — 3 for 4.5 < z < 6. As a consequence, the UV 



metagalactic background shortward of the Lyman edge may 
have undergone a qualitative change in its spectral shape in 
the interval 3.5 < z < 4.5. QSOs produce too few hydro- 
gen ionizing photons to reionize the IGM any earlier than 
z T i w 4, even assuming a hard spectral index (as = 0.5). 
For a moderate spectral index (as = 1.8), QSOs alone are 
unable to ionize the helium in the IGM any earlier than 
z ~ 3, but helium reionization may occur as early as z ~ 5 
for QSOs with a hard spectral index (as = 0.5) (unless the 
clumping factor is large). If QSOs do not ionize He n until 
z ?s 3, then a second qualitative change will result in the 
spectral shape of the UV background. 

Both PLE and PDE models are found to provide accept- 
able fits to the QSO counts over the redshift range 3 < z < 6. 
Comparison with PLE models for z < 2, however, suggests 
substantial density evolution in the QSO luminosity func- 
tion between z > 3 and z ^,2. Predicted break magnitudes 
for the models at z > 3 suggest that the breaks would be 
clearly detectable in the PDE case for a survey one magni- 
tude deeper than the SDSS QSO survey and three magni- 
tudes deeper for the PLE case. 

Consistency with the WMAP 2a lower limit of z T i > 11 
(Kogut et al. ) suggests that the sources which reionized 
the IGM must have had hard spectra (as « 0.5), like hard 
spectra QSOs or AGN, or Pop III stars, unless reionization 
occurred very early (z » 6), in which case the ionization 
state of the IGM at z < 6 would not be expected to yield 
any constraints on the nature of the reionization sources. 
An assessment of the contribution of ionizing photons from 
QSOs based on recent surveys shows that they fall short of 
the minimal requirement for reionizing the IGM by a factor 
of several. This suggests that either: 1. galaxies were much 
more prodigious sources of ionizing photons at z > 6 than 
at z < 6, either because of a higher star formation rate or a 
higher escape fraction of ionizing photons, 2. star formation 
in galaxies at z > 6 was dominated by the production of 
Pop III stars, producing a hard spectrum, 3. at least a small 
fraction of galaxies harboured a low luminosity AGN at z > 
6, so that they were hard spectrum sources at these epochs, 
or 4. reionization was produced by an as yet unidentified 
population of sources (such as miniquasars). 

Another possibility is that numerical simulations of the 
structure of the IGM at high redshifts are seriously in er- 
ror. The high precision agreement with measurements of the 
Lya forest in spectra at moderate redshifts (z < 4) (Meiksin, 
Bryan & Machacek 2001) does not preclude the possibil- 
ity that energetic events such as galactic winds or QSOs 
at higher redshifts substantially disturbed the gas at ear- 
lier times, so that current simulations do not provide valid 
constraints on the source emissivity required to match the 
mean Lya transmitted flux measurements at z > 5. Cur- 
rently there is no direct evidence that this may be the case. 
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